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Abstract 



Using a classical analjd;ical criterion (that of curvature) and numerical 
results (Poincare sections and spectral statistics), a transition order-chaos- 
order in the roto-vibrational model of atomic nuclei has been shown. Nu- 
merical calculations were performed for some deformed nuclei. 



PACS numbers: 21.10.-k, 05.45.+b, 03.65 Sq 
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1. Introduction 

In the last few years many authors have shown great interest in the so- 
called "quantum chaos": the properties of quantal systems which are chaotic 
in the semiclassical hmit — > [1,2,3]. 

In atomic nuclei the coexistence of regular, chaotic and collective states 
makes the problem quite intricate. In order to disentangle it, many models 
have been used (see, for example [4,5,6] and references quoted therein). The 
aim of this paper is to study the transition from ordered to chaotic states in 
the roto-vibrational model of atomic nuclei. This model was introduced by 
Bohr and Mottelson [7] and discussed in great detail in [8,9,10]. 

2. The Model 

Since the model has been amply described in [8], in this section we limit 

ourselves to reporting only a few basic formulae. 
The hamiltonian is: 

H = E^if, + Ej-ot: (1) 

where: 

E^ib = ^B{al + 2al) + V{ao, 02), (2) 
1 3 ^ 

Erot = 9 m ^k'^k{(^0, «2), (3) 
^ fe=l 

with: 

V{ao, as) = E Cmn(a^ + 2al)'''a^{6al - a^)". (4) 

m,n 
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The parameters ao and a2 are connected to the deformation (3 and asymmetry 
7 by the standard relations [8]: 

ao — Pcos'j, 02 = -7= sin 7. (5) 
v2 

In terms of the new variables the components of the moment of inertia are 
[7]: 

27r 

J, = 4S/32sin2(7-yfc). (6) 
If the nucleus has an axially symmetric deformation: 

UJ 

uji = uj2 = a;3 = 0, (7) 

and if we take only the first terms of (4), the hamiltonian (1) can be written: 
H = \B{al + 2al) + V{ao, 02) + \BiJ^{?>al + 2a^), (8) 

where: 

y(ao, 02) = ]f2{al + 2a2) + ^CsaoiGaj - al) + + 2«2)' + ^0, (9) 

with Vq a parameter chosen to have the minimum of the potential equal to 
zero. As discussed in [8], the presence of bound states in atomic nuclei leads 
to a value of C4 > 0, whereas for C3 a positive value corresponds to a prolate 
shape, a negative value to an oblate shape. Similarly C2 may also be either 
positive or negative. 

The shape of the nuclear potential V{ao, 02) is a function of C2 and x — 
Cl/{C2Ci). For C2 > 0, and < x < 56/9 the nucleus is spherical; for 
56/9 < X < 7 the nucleus is spherical in the ground state (g.s.) and deformed 
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in the excited states (e.s.); for x > 7 it is deformed in the g.s. and spherical 
in the e.s; for C2 < it is deformed in the g.s. and 7-unstable in the e.s. 
Near the equihbrium (do = /3o, 02 = 0), y(ao, 02) can be written: 

V(ao, 02) = Y(ao - Pof + C^al (10) 

In order to calculate the equilibrium configuration the following system 
must be solved: 

dV 4 

= /3o(C2 - 3y ^C^3/9o + -^C^al) = 0, (11) 

^1/30,0 = C2 - 6y-C3/3o + yC^4/3o = ^l/3o,0 = Co, (12) 

= + ^/l^^^" + f = l^''^-^ = ^^^^ 

In roto-vibrational nuclei the parameters and 6*2 are simply connected 
to some experimental quantities: 

Ep = fiup = = huj^ = h\l (14) 



and: 



^°-3^' ^^"3^' ^^^^ 
where e = fi^ j Jo is taken from the experimental energy of the first rotational 
state and /3o is the equilibrium deformation of the nucleus, which may be 
obtained from the 5(^2, 2+^0+) value [8]. Finally: 



1 

18e/?o' 



^2 = -T^(3£;|-^^), 
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As a simple application of the above procedure, the nuclear potentials 
have been calculated for ^^*Dy and ^^^Er (sec Figures 1,2). The numerical 
values of the parameters /Sq, e, Ep, E-y, taken from the experimental data, 
are shown in Table 1. 

3. The onset of chaos by curvature criterion 

As is well known, the transition order-chaos in systems with two degrees 
of freedom may be studied by the curvature criterion [11,12]. It is however 
important to point out that in general the curvature criterion guarantees 
only a local instability and should therefore be combined with the Poincare 
sections [14]. For a fuller discussion of this point see [13]. As mentioned in 
the previous section, the total hamiltonian can be written: 

H =^B{al + 2dl) + W{ao,a2), (17) 

where: 

W{ao, as) = V{ao, as) + ^5cu^(3a^ + 2al) (18) 
is the effective potential. As a function of Cs, C3, C4, (18) becomes: 

W{a^, as) = ^(:72(a^+2a^)+y^C3ao(6a2-ag)+lc4(ag+2a2) Vyo+^5c^'(3a2+2a^). 

(19) 

Owing to the symmetry properties of the effective potential W ^ our study 
may be restricted to the case H^(ao, as = 0) [8]. To apply the above criterion 
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to our system, the sign of the curvature K can be studied by solving the 
equation: 

d'^W , 12 „ o „ / 2 

whose discriminant A is given by: 



K{ao) = ^(ao, 02 = 0) = ^C^al - ^sJ^C^ao + (C2 + Wuj'') = 0, (20) 



A = |c2C.(x-y-14^). (21) 

If A < the curvature K is always positive and the motion is regular; if 
A > there is a region of negative curvature and the motion may be chaotic. 

For C2 > 0, < X < 14/3 (spherical nuclei) the curvature is positive 
and therefore the motion is regular for all uj. For % > 14/3 (spherical and 
deformed nuclei) and: 



the curvature is negative and chaotic motion may appear. 

For C2 < (7-unstable nuclei) there is a region with negative curvature 
for < a; < ujci where: 



is the critical frequency of the system. As a function of e and /3o the critical 
energy can be written: 



14 3e^ 

= Y (X - y ) ^4 (24) 



and the critical angular momentum Ic is: 



It is perhaps noteworthy that the shape of the effective potential W{ao, 0) 
changes drastically as a function of lu (see Fig. 3). If a; increases there is a 
transition from chaos to order: the region of chaotic motion decreases and 
becomes zero for cu > lOc- For a fixed value of the parameter x, the region of 
chaotic motion is limited, on the line A'(ao) = 0, by the two branches defined 

Table 3 shows, for some nuclei, the energies -Ec/taos ^"^^ ^thaos^ which limit 
the chaotic region for a; = 0, and also the critical angular momentum /c- 

4. Numerical study of the order— chaos transition 

As mentioned in the previous section, the curvature criterion is able to 
characterize the local behaviour of the system (for example the local insta- 
bility) and may give only a rough signature of the global properties (e.g. the 
global instability) [13]. As is well known, a very useful tool for the study 
of global properties is provided by the Poincarc sections [14]. With this 
aim the classical trajectories have been calculated by a fourth order Runge- 
Kutta method. In order to avoid numerical errors connected to the use of 
finite temporal intervals, a first-order interpolation has been used [15] . 

The Hamilton equations of the systems are: 

do = -Bpo, 

Po = -Csao - 2^C3(3a^ - 3a^) - ^C^ao{al + 2al) - Wu'^ao, (27) 
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P2 = -2C2a2 - 12\ — 




6*30002 - -C402(oo + 202) - 2BuP'a2, 



where po and p2 are the conjugate momenta: 



Po = -Boo, P2 = 2-602. 



(28) 



Figure 4 shows the Poincare sections for ^^^Gd at the energy 5.5 MeV and 
for different values of rotational frequency. The figure clearly shows a chaos- 
order transition as the frequency u increases. In Figure 5 for ^^^Er the 
Poincare sections are shown for different values of the energy and rotational 
frequency a; = 0. As can be seen, there is a chaos-order transition, albeit 
not so sharp as in the previous case. 

It is well known that the fluctuation properties of quantal systems with 
underlying classical chaotic behaviour and time-reversal symmetry are in 
agreement with the predictions of the Gaussian Orthogonal Ensemble (GOE), 
and that quantum analogs of classically integrable systems display the char- 
acteristics of Poisson statistics [4,6]. 

For deformed nuclei, like those of the rare-earth region, it is not easy to 
obtain the correct energy levels, because the potential energy is an asym- 
metric triple well. To avoid the problems related to the tunneling effects, 
a numerical method based on the formulation of quantum mechanics using 
the cuclidcan path integral [16] should be used. The energy spectrum is 
then mapped into one with quasi-uniform level density by means of the local 
unfolding procedure described in detail in reference [17]. 

In Figure 6 the spectral statistics P{s) and A3 [18,19] are plotted for 
^^^Er. These statistics confirm the classical results: for energies above the 
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saddle energy, about 4 MeV, there is prevalently chaotic behaviour; for higher 
energies, about 15 MeV, there is mixed behaviour with a predominance of 
regular classical trajectories. The non-universal behaviour of A.s{L) for large 
values of L, not predicted by GOE, has been explained by Berry [20] using 

the semiclassical quantization. 

5. Conclusions 

Wc have shown by combining analytical results (the curvature criterion) 
with numerical ones (the Poincare sections and spectral statistics) that, in 
the roto-vibrational model of atomic nuclei, an order-chaos-order transi- 
tion occurs as a function of the energy. Our results are in good agreement 
with those of [12], but, unlike the authors of [12], we have also calculated 
the energy ranges of the chaotic regions (see Table 3) and the chaos-order 
transition which occurs as a function of the rotational frequency lu. 

The authors are greatly indebted to M. Rosa-Clot and S. Taddei for the 
provided numerical data of energy levels. 
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TABLE CAPTIONS 



Table 1: Numerical values of the parameters (3o, e, Ep, E^, taken from the 
experimental data for some even-even nuclei. The parameters are defined in 
the text. 

Table 2: Numerical values of the parameters C2, C3, C4, Vq defining the 
effective potential W. 

Table 3: Numerical values of the energies which delimit the chaotic region 
for a; = and the critical angular momentum. 



13 



FIGURE CAPTIONS 



Figure 1: Nuclear potential for Dy. 
Figure 2: Nuclear potential for ^^^Er. 

Figure 3: Effective nuclear potential for ^^'^Dy as a function of the rota- 
tional frequency co; from left to right: hu; — MeV, fiu; — 0.5 MeV, ku; — 1 
MeV. 

Figure 4: The Poincarc sections for ^^^Gd at the energy 5.5 MeV and for 
different values of rotational frequency; from the top: hcu — MeV, hcu — 0.5 
MeV, hu = 1 MeV. 

Figure 5: The Poincare sections for ^^^Er at the rotational frequency 
cl; = and for different values of the energy: (a) E = 1 MeV, (h) E = 6 
MeV, (c) E = 9 MeV, (d) E = 12 MeV. 

Figure 6: Spectral statistics P{s) and A^lL) for ^^^Er at the rotational 
frequency a; = for different energy regions: 2 < E < 6 MeV (below) and 
for 13 < £; < 17 MeV (above). The sohd hne is the GOE statistic curve and 
the dashed line is the Poisson one. 
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Nucleus 


Po 


e (MeV) 


Ef3 (MeV) 


(MeV) 




0.47 


0.022 


1.80 


0.96 




0.41 


0.021 


1.76 


0.73 


166^^ 


0.33 


0.023 


1.14 


0.75 




0.23 


0.015 


0.63 


0.76 


238 JJ J. 


0.28 


0.014 


0.99 


1.05 


240p^ 


0.28 


0.013 


0.86 


0.92 



Table 1 
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Nucleus 


C2 (MeV) 


C3 (MeV) 


C4 (MeV) 


Vo (MeV) 




-100.64 


37.57 


668.01 


5.56 




-138.41 


34.97 


1155.67 


5.69 


166 jjjj. 


-73.77 


64.96 


1138.52 


1.89 




-43.56 


229.67 


2260.62 


0.44 


238 JJ J. 


-94.54 


334.20 


3276.70 


1.56 


240p^ 


-73.64 


281.53 


2664.58 


1.20 



Table 2 
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Nucleus 


Kkaos (MeV) 


Khaos (MeV) 


icih) 




2.36 


5.36 


31.86 




2.51 


5.59 


33.31 


166 


0.83 


1.89 


18.79 




0.15 


0.44 


13.28 


238 JJ J. 


0.61 


1.55 


23.88 


240p^ 


0.46 


1.19 


21.97 



Table 3 
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